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SECTION  I 


INTRODUCTION 

Vibrations  caused  by  acoustic  pressure  can  frequently  disturb  the  operat¬ 
ing  conditions  of  various  instruments  and  systems,  and  sonic  fatigue  failures 
which  occurred  in  aircraft  structural  components  cause  large  maintenance  and 
inspection  burdens  for  the  Air  Force.  The  development  of  sonic  fatigue  data 
and  design  techniques  were  initiated  to  prevent  sonic  fatigue  failures. 

Design  methods  and  design  criteria  for  many  types  of  aircraft  structures  have 
been  developed  under  Air  Force  sponsorship  and  by  the  industry  in  the  past 
twenty  years.  Reference  1  has  a  complete  list  of  the  reports  describing 
these  efforts.  This  research  led  to  sonic  fatigue  design  criteria  and  design 
charts  which  are  widely  used  during  the  design  of  an  aircraft.  Although  cur¬ 
rent  analytical  sonic  fatigue  design  methods  are  essentially  based  on  small 
deflection  or  linear  structural  theory  (Reference  1,  page  209),  many  docu¬ 
mented  tests  on  various  aircraft  panels  have  indicated  that  high  noise  levels 
produce  nonlinear  behavior  with  large  amplitudes  in  such  structural  panels. 

For  example,  Fitch  et  al ,  (Reference  2),  van  der  Heyde  and  Smith  (Reference  3), 
Jacobs  and  Lagerquist  (Reference  4),  and  Jacobson  (References  5  and  6)  have 
repeatedly  reported  that  a  poor  comparison  exists  between  the  measured  and 
computed  root  mean  square  (RMS)  displacements  and/or  RMS  stresses.  They  all 
observed  that  the  test  panels  responded  with  large  deflections  at  high  sound 
pressure  levels,  whereas  the  analytical  responses  were  based  on  linear  small 
deflection  theory.  The  neglect  of  such  large  deflection  geometrical  non¬ 
linearity  in  analysis  and  design  formulations  has  been  identified  as  one  of 
the  major  causes  for  disagreement  between  experimental  data  and  analytical 
results.  The  evidence  of  those  researchers  was  summarized  in  Reference  7, 
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where  a  comprehensive  review  of  existing  analytical  methods  on  random  exci¬ 
tations  of  nonlinear  systems  was  also  given. 

Because  there  are  no  reliable  analysis  methods  available  for  predict¬ 
ing  the  nonlinear  stress-sound  pressure  relation,  costly  and  time-consuming 
full-scale  fatigue  tests  of  aircraft  structures  and  components  are  frequently 
conducted.  The  objectives  of  the  present  work  are:  1)  to  gain  a  better 
understanding  of  the  random  response  of  nonlinear  panels,  2)  to  incorporate 
large  deflection  geometrical  nonlinearity  into  analysis  methods  for  determin¬ 
ing  structural  response  to  high  intensity  noise,  and  3)  to  provide  analytical 
background  material  for  formulation  of  improved  sonic  fatigue  design  proce¬ 
dures  that  would  result  in  better  and  less  costly  designs  without  sacrificing 
safety. 

The  Karman-Herrmann  large  deflection  equations  for  rectangular  plates 
(Reference  8)  are  employed  in  this  development.  Using  a  single-mode  Galerkin's 
approximation,  the  dynamic  equations  reduce  to  a  nonlinear  differential  equa¬ 
tion  with  time  as  the  independent  variable.  The  method  of  equivalent  lineari¬ 
zation  is  then  applied  to  reduce  the  nonlinear  equation  to  an  equivalent 
linear  one  (References  12-14).  Mean-square  displacements,  mean-square  stresses, 
and  nonlinear  frequencies  at  various  acoustic  loadings  are  obtained  for  rec¬ 
tangular  panels  of  different  aspect  ratios  and  damping  factors.  Both  simply 
supported  and  clamped  boundary  conditions  with  immovable  and  movable  inplane 
edges  are  considered.  The  results  are  presented  in  graphical  form.  Compari¬ 
sons  with  experimental  results  are  also  presented. 
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SECTION  II 


FORMULATION  AND  SOLUTION  PROCEDURE 
1.  Governing  Equations 

Assuming  that  the  effect  of  both  the  in  plane  and  rotatory  inertia  forces 
can  be  neglected,  the  dynamic  equations  of  a  rectangular  isotropic  plate 
undergoing  moderately  large  deflections  are  (References  8  and  9): 

L(w,F)  =  Dv4w  *  phw,tt  -  h  <F,yy  w,xx  ♦  F,xx  w,yy 

"2F,Xy  W,Xy  )  -  p(t)  =0  (1) 

v4F  =  E  (w2,xy  -  w,xx  )  (2) 

where  w  is  the  transverse  deflection  of  the  plate,  h  is  the  panel  thickness, 
p  is  the  mass  density  of  the  panel  material,  D  =  Eh  /12(l-o  )  is  the  flexural 
rigidity,  E  is  Young's  modulus,  j  is  Poisson's  ratio,  p(t)  is  the  exciting 
pressure,  and  a  comma  preceding  a  subscript(s)  indicates  partial  differentiation  (s). 
The  stress  function  F  is  defined  by 

°x  =  F\yy 

oy  =  F,xx.  (3) 

Txy  ”  ~F’xy 

where  o  ,  o  ,and  t  are  membrane  stresses, 
x  y  xy 

a.  Simply  Supported  Panels 

For  a  rectangular  plate  simply  supported  along  all  four  edges  as  shown  in 
Figure  1,  Chu  and  Herrmann  (reference  8),  an-'  Lin  (reference  ln)  have  consi¬ 
dered  that,  if  the  f’l-iarentrl  nedr  is  predominant,  the  notion  of  the  panel  can 
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be  represented  adequately  as 


w  =  q(t)  h  cos  (it  x/a)  cos  (w  y/b)  (4) 

where  q(t)  is  a  function  of  time  only.  The  maximum  value  of  q(t)  coincides 

with  the  maximum  deflection  wmau  divided  by  panel  thickness  h.  The  expression 

max 

w  satisfies  the  boundary  conditions  for  simple  supports. 


w  =  w’xx  +  v  w,yy  =  0.  on  x  =  ±  a/2 


w  =  W’yy  +  v  w,xx  =  0*  on  y  =  ±  b/2 


(5) 


Substituting  the  expression  for  w  in  Eq.  (2)  and  solving  for  a 
particular  solution  Fp  yields 

Fp  =  -(3})  q2h2Er2  (cos  ^  ^  cos  (6) 

where  r  =  a/b.  The  complementary  solution  to  Eq.  (2)  is  taken  in  the  form 

Fc  ■  \  4  +  fy  4  -  "xy 

where  the  constants  FT,  FT,  and  FT  contribute  to  the  membrane  stresses 


(7) 


V  °y  and  Txy  and  are  t0  be  determined  from  the  inplane  boundary,  irrmovable 
or  movable,  conditions. 

For  the  immovable  edges  case,  the  conditions  of  zero  inplane  normal 
displacement  at  all  four  edges  are  satisfied  in  an  averaged  manner  as 


vF’xx}  -  *  w 


’X 


]  dxdy 


,on  x 


+  a/2 


(8) 

on  y  *  +b/2 


where  u  and  v  are  inplane  displacements. 

For  the  movable  edges  case,  the  edges  are  free  to  move  as  a  rigid 
body  with  the  average  inplane  stress  equal  to  zero.  The  inplane  boundary 
conditions  are 


F-  *  -fln 

"* '  %2  W  - 0 


on  x  =  +a/2 


u  =  constant 


F,  =0 
n/l 

N-  -  hi  F,  uw  dx  =  0 
J-a/2 


yy 


on  y  =  ±b/2 


v  =  constant 


(9) 


where  N  and  N  are  membrane  stress  resultants  per  unit  length  in  plate. 
a  y 

By  making  use  of  these  inplane  edge  boundary  conditions,  Eqs.  (8)  and 
(9),  it  easily  can  be  shown  that  for  the  immovable  edges 
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2.2,  2 


N„  ’  q  2  >  (1  *  vr‘) 

x  8a'(l-v  ) 


2 1  2  2 
_  q  h  En 


y  8a2(l-v2) 


( r  +  v) 


00) 


*xy 


and  for  the  movable  edges 


"*  ■  "y  ■  Nxy  ■  0 


(11) 


the  complete  stress  function  is  then  given  by  F  =  Fp  +  Fc> 


With  the  assumed  w  given  by  Eq.  (4)  and  stress  function  given  by 
Eqs.  (6)  and  (7),  Eq.  (1)  is  satisfied  by  applying  Galerkin's  method 


a/2  ,b/2 


J  J 


L  (w,F)  w  dxdy  =  0 


(12) 


-a/2  •'-b/2 

from  which  yields  the  modal  equation  of  the  form 

and 


(13) 


2.  ,2  D 
X°  ^  ’ 


,2  4  hJ  \2 

X- =  it  (1  +-jr) 
r 


m  =  ir2ph2/16 


(14) 


e  *6p+  8C 


1 


with 


(1  ♦  r4)  (1-v2) 
4r 


B 


c 


* 


[1  ♦  vr2  +  r2(r2  ♦  v)] 
Zrq 


where  w  is  linear  radian  frequency,  mis  mass  coefficient,  and  B  is  non- 

o 

linearity  coefficient.  The  linear  frequency  XQ,  nonlinearity  coeffi¬ 
cients  b!  and  2*,  and  aspect  ratio  r  are  all  nondimensional  parameters. 

r  C 

b.  Clamped  Panels 

Yamaki  (Reference  11)  considered  the  predominant  mode 


w  =  al|)h  (1  +  cos  2£X)  (1  +  cos  (IS) 

which  satisfies  the  clamped  support  conditions 

w  =  w,x  =  0  ,  on  x  =  ±  a/2  (16) 

w  =  w,y  =  0,  on  y  =  ±  b/2 

By  introducing  Eq.  (15)  in  Eq.  (2)  and  solving  it,  the  particular 
stress  function  is 


2.  2_  2 

q  h  Er 


cos  +  IF  C0S 


4ttx 

a 


+ 


(1+  r2)2 


(4+r2)2 


cos 


4ux 


cos^* 


cos 


2nx 

a 


cos 


4-ny 


(17) 
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The  complementary  stress  function  Is  assumed  as  the  form  appearing 
in  Eq.  (7).  Upon  enforcing  the  Inplane  edge  conditions,  Eqs.  (8)  and 
(9),  it  can  be  shown  that  for  the  immovable  edges 


(1  +  vr  ) 


*  3g2h2E*2 
32  aZ(l-v2) 

•  3#*2  r  (r2+v) 
32  a^(l-v2) 


(18) 


Nxy  =0 


for  the  movable  edges 

Nx  =  Ny  =  Nxy  =  0  (19) 

the  complete  stress  function  is  given  by  F  =  Fp  +  Fc>  Introducing 
these  expressions  for  w  and  F  in  Eq.  (1)  and  applying  Galerkin's  proce¬ 
dure  yields  the  modal  equation 

q  +  J  q  +  eq3  =  ^  (13) 


where 


2  _  .2  0 
O)  -  A  - 1 

0  0  Phb^ 


(3  +  2r2  +  3r4) 
9r 


m  =  9  ph  /16 


(20) 


6  =8  +  3  =  (e*  +  3*)  — =t r 
p  C  P  c 


8 


and 


2 


4 

TT 


(1-v2) 


(1 


+  r2)2 


2(4+r2)2  2(l+4r2)2  1 

g*  =  —  [  1  +  vr2  +  r2  (r2  +  v)] 

c  2r4 


(21) 


Equation  (13)  represents  the  undamped,  large  amplitude  vibration  of  a  rec¬ 
tangular  panel  with  simply-supported  and  clamped  edaes. 

The  methods  commonly  used  for  determining  the  damping  coefficient  are 
the  bandwidth  method  in  which  half-power  widths  are  measured  at  modal  reso¬ 
nances  and  the  decay  rate  method  in  which  the  logarithmic  decrement  of  decay¬ 
ing  modal  response  traces  is  measured.  The  values  of  damping  ratio 
?  C/C(J  range  from  0.005  to  0.05  for  the  common  type  of  panel  construc¬ 
tion  used  in  aircraft  structures.  Once  the  damping  coefficient  is  determined 
from  experiment;  or  from  existing  data  of  similar  construction,  the  modal 
equation,  Eq.  (13),  now  reads 


V  *  2w„q  ♦  »02q  *  Bq3  .  EM  (k) 

The  method  of  equivalent  linearization  iz  employed  to  determine  an  approxi¬ 
mate  RMS  displacement  from  Eq.  (22). 

2.  Method  of  Equivalent  Linearization 

The  basic  idea  of  the  equivalent  linearization  (Reference  12-14)  is  to 
replace  the  original  nonlinear  equation,  Eq.  (22),  with  an  equation  of  the 
form 


9 


(23) 


w 


q  +  2  c«0q  +  ft2q  +  err  (q)  = 

where  n  is  an  equivalent  linear  or  nonlinear  frequency,  and  err  is  the 
error  of  linearization.  An  equivalent  linear  equation  is  obtained  by 
omitting  this  error  term,  then  Eq.  (23)  is  linear  and  it  can  be  readily 
solved.  The  error  of  linearization  is 

err  =  (u)Q2  -  ft2)  q  +  Bq3  (24) 

which  is  the  difference  between  Eq.  (22)  and  Eq.  (23).  The  equivalent  linear 
frequency  J2  is  chosen  in  such  a  way  as  to  make  the  error  of  linearization 
term,  err  (q),  as  small  as  possible.  To  this  end  the  mean-square  error 

err2  is  minimized  ,  that  is  3(err^)  _  n  (25) 

HnZ) 

If  the  acoustic  pressure  excitation  p(t)  is  stationary  Gaussian  and 

ergodic,  then  the  response  q  computed  from  the  linearized  equation, 

Eq.  (23),  must  also  be  Gaussian.  Substituting  Eq.  (24)  into  Eq.  (25) 

yields  (References  10  and  12) 

2  2  ~T 

ft  =  «o  +  3  Bq  (26) 

~7 

where  q  is  the  maximum  mean-square  deflection  of  the  panel.  Dividing 
both  sides  of  Eq.  (26)  by  D/phb^  yields 

2  2  *  ~7 

A  *  *0  +  33  q‘  (27) 

2 

where  A  is  a  nondimensional  equivalent  linear  or  nonlinear  frequency 


parameter. 
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An  approximate  solution  of  Eq.  (23)  is  obtained  by  dropping  the 
error  term,  the  mean-square  response  of  modal  amplitude  is 

S(aj)l  H<^2  da,  (28) 

where  S  (w)  is  the  spectral  density  function  of  the  excitation  pressure 
p(t),  and  the  frequency  response  function  H(oj)  is  given  by 

HU)  •  J-  -  <29> 

m(ft  -  u  +  2 i ) 

For  lightly  damped  (c  f_0.05)  structures,  the  response  curves  will  be 
highly  peaked  at  n.  The  integration  of  Eq.  (28)  can  be  greatly  simplified 
if  the  forcing  spectral  density  function  S(w)  can  be  considered  to  be 
constant  in  the  frequency  band  surrounding  the  nonlinear  resonance  peak 
fl,  so  that 

7*  -4^  (30) 

q  4nrcu>  $T 
0 

In  practice,  the  spectral  density  function  is  generally  given  in  terms 
of  the  frequency  f  in  Hertz.  To  convert  the  previous  result  one  must 
substitute 

n  =  2  nf 

and  S(n)  =  S(f)/2»  (31) 


into  Eq.  (30),  the  mean-square  peak  deflection  is  simply 


32  Sf 

- - —  t  for  simply  supported  panels 

T  »  a  A 

q  m  ^  0 

32  Sf 

81cX  a* 
o 


for  clamped  panels 


(32) 


The  pressure  spectral  density  function  S(f)/2it  has  the  units 
2  2 

(Pa)  /Hz  or  (psi)  /Hz,  and  Sf  is  a  nondimensional  forcing  exci¬ 
tation  spectral  density  parameter  defined  as 


c  _  s(f) 

f'  .WF 


(33) 


The  linear  frequency  parameters  XQ  in  Eqs.  (32)  are  given  in 

Eq.  (14)  and  Eq.(20)  for  simply  supported  and  clamped  panels, 

respectively,  and  the  equivalent  linear  frequency  parameters 
2 

x  can  be  determined  through  Eq.  (27). 

3.  Solution  Procedure 


The  mean-square  response  q  in  Eq.  (30)  (Or  Eq.  (32))  is 
determined  at  the  equivalent  linear  frequency  0  (or  x)  which  is  in 
turn  related  to  q  through  Eq.  (26)  (or  Eq.  (27)).  To  determine  the 
mean-square  deflection,  an  iterative  procedure  is  introduced.  One  can 
estimate  the  initial  mean-square  deflection  qQ  using  the  linear  frequency 
through  Eq.  (30)  as 


2  3 

4m  s os 

o 


(34) 


This  initial  estimate  of  qQ  is  simply  the  mean-square  response  based 

~7 

on  linear  theory.  This  initial  estimate  of  qQ  can  now  be  used  to  obtain 

2  2  ~ 2" 

refined  estimate  of  Jlj  through  Eq.  (26),  n,  =  +  36  qQ  ,  then 

q^  is  obtained  through  Eq.  (30)  as 


q?  =  — 

4m  CWqOj 


(35) 
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As  the  iterative  process  converges,  the  relation 


_ *  S(nn) 

2  2 
4m  o  _ 
o  n 


(36) 


becomes  satisfied.  In  the  numerical  results  presented  in  the  following 
section,  convergence  is  considered  achieved  whenever  the  difference  of  the 
RMS  displacements  satisfied  the  relation 

<  10'3  (37) 


4.  Stress  and  Strain  Response 

Once  the  RMS  nodal  displacement  is  determined,  the  bending  stresses 
on  the  surface  of  the  panel  can  be  determined  from 


°xb  =  'fT  <w’xx  +  v  w’yy} 

(38) 

°yb  =  (w’yy  +  v  w’xx) 


The  corresponding  strains  are  given  by 

e  -  -  few 
exb  2  w’xx 

h 

eyb  "  7  w’yy 


(39) 


The  membrane  stresses  in  the  panel  are  obtained  from  Eqs.  (3) 
and  the  corresponding  stress  function  F.  Membrane  strains  are  given  by 


exm  =  E  ^F’yy  "  vF’xx^ 


ym 


*  ^  (F» 


xx  "  vF’yy^ 


(40) 
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a.  Simply  Supported  Panels 

From  Eqs.  (3)  and  (38),  and  using  Eqs.  (4),  (6),  (7)  and  (10),  the 
expressions  for  the  nondimensional  stresses  on  the  surface  of  the  panel 
with  immovable  edges  are  given  by 


(  0  t  +  0  )  - ry 

v  xb  xm'  £^2 


h('\lry  (  L2  +  v)  cos  IT  cos  ^b1 

'  r 


+  (-4  cos  &-)  q2  +[  4"  v-r-k 


8r  (1 -v2) 


a  bc  h2 

■  <>  *  >> 


l-~  (1  +  -2 

2(l-v2)  r‘ 


(1  +  -2)  cos  jf- cos  j^-]  q 


2 (r2+v)  ,  2 


Hf  cos  -2J2.  )  q2  +[  q 

8  a  8r  (1-v2) 


The  tensile  strains,  Eqs.  (39)  and  (40),  on  the  surface  of  the  panel 
are  then  given  by 


_r  tt2  / h \2 

ex  ~  (b}  C0S 


a  'os  bL]  q 


+  \  cos  nF  ■ v  cos  in ) q 

r 

V yr* ' Ar*  * v>  j  <i2 

8  b  rZ(l-v!) 


2itx  ,  „2 
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(42) 


(42) 


eyS(  cos  ?  cos  1  q 

+  r<F)2  (cos  H5- "  -2  cos  ^  ^2 

r 

♦  AV  r 2  *  V'l(1./>  ur2)1  q2 

W  717-7) 

For  movable  inplane  edges,  the  last  term  in  Eqs.  (41)  and  (42)  vanishes, 
b.  Clamped  Panels 

Similarly,  from  Eqs.  (3)  and  (38)  and  using  Eqs.  (7),  (15),  (17),  and 
(18),  the  expressions  for  the  nondimensional  tensile  stresses  on  the  sur¬ 
face  of  a  clamped  panel  with  immovable  edges  are 


°  b  2  1 

— =  - -  [  cos  a  x(l  +  cos  B  y)  +  v(l+cos  A  x)  cos  B  y]  q 

Eh  2(1 -v  )  r 

2  2 

+  - -p-  (V  COS  B  y  +  - Y  cos  A  x  cos  B  y  +  -^r  cos  2  B  y 

8  r  (1+rV  4r4 

+ - K-tj  cos  2  A  x  cos  B  y  +  - - - cos  A  x  cos  2  B  y]q2 

(4  +  r  r  (1  +  4r  ) 


+  [^1+-vf>i  q2 

32r^(l -v  ) 


(43) 


o.  b 


-^5—  =  - 2— __  [(1  +  cos  a  x)  cos  B  y  +  \  cos  A  x  (  1+cos  B  y)]  q 

ElT  2(1V)  r 


v2  2 

+  s  [cos  A  x  +  k  cos  2  A  x  +  - cos  A  x  cos  B  y 

8  (1  +  ry 

4  *  r 

+  - 5-tj-  cos  2  A  x  Cos  B  y  +  - 5—5-  cos  A  x  cos  2  B  y]  q‘ 

(4+r^r  (l+4r^r 

♦  c3'-Y2ty|  1  q2 

32r‘(l  V) 
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where  A  =  2  n/a  and  B  =  2  n/b.  The  corresponding  '.ensile  strains  are 
given  by 

ex  =  (5-)^  cos  A  x  (  1  +  cos  B  x)]  q 

2r 

2  2 

.  n  r  ,h>2  fl  D  2 

+  '  a  "%)  [“A  cos  B  y  +  - 5-9-  cos  A  x  cos  B  y 

r  (1  +  r2)2 

+  7T  cos  2  B  *  +  T — hv2  cos  2  A  x  cos  B  y  +  - ^-*-x  cos  A  x  cos  2  B  y] 

4r  (4+r*r  (l+4r2)2 

_  _v  [cos  fl  x  +  ^  cos  2  A  x  +  - - — cos  A  x‘  cos  B  y 

r2  (1  +  r2)2 


4  1  r 

+  - 272  cos  2  A  X  cos  B  y  +  - ± - «-*-  cos  A  x  cos  2B  y]  q 

(4  +  r£)£  (1  +  4r2)2 


4.  rll2  /h^2  1  +  vr2-v(r2+v)  n  2 

32  ?  t -  1  q 

r  (1-v  ) 


ey  “  tyr  (jf)^  (  1  +  cos  A  x)  cos  B  y]  q 


t2  ,h,2 


(£)2  jcos  A 


x  +  h  cos2A  x  +  - —x  b  cos  A  x  cos  B  y 

(1  +  r2)2 


4  1 

+  7  2v2  cos  2  ^  x  cos  B  y  + - -9-9-  cos  A  x  cos  2  B  y 

(r  +  r‘r  (1  +  4r2)2 


t“4  cos  B  y  +  - cos  A  x  cos  B  y  +  cos  2  B  y 

I  -I  .4  ^ 


(l+r‘) 


+  TT^W  cos  2  A  x  cos  B  y  +  - -  x  9  cos  A  x  cos  2B  y]  l  qZ 

H+r'T  ( l+4r2)2  / 


4.  r 3tt2  /h^2  r2+v-v(l+vr2),  2 

1  32  {b}  —JZ  2 x 


'  (1 -V  ) 


For  movable  edges,  the  last  term  in  Eqs.  (43)  and  (44)  vanishes. 

Examining  Eqs.  (41)  -  (44),  a  general  expression  is  obtained  for  the 
stress  (or  strain)  at  any  point  in  the  structure  as 

0  =  0^  +  C2q2  (45) 

where  and  C2  are  constants.  The  expressions  for  and  can  be  found 

from  Eqs.  (41)  to  (44).  The  constants  can  be  determined  from  material  pro¬ 
perties,  dimensions  of  the  panel,  and  the  location  and  direction  at  which 
the  stress  is  to  be  measured.  The  mean-square  stress  (or  strain)  is  then 
related  to  the  mean-square  modal  amplitude  in  a  general  expression  as 

a2  -  C2  q7  +  3  C2  (?)2  (46) 

2 

Once  the  mean-square  deflection  q  is  determined,  Eqs.  (36)  and  (37),  the 
mean-square  stress  (or  strain)  can  then  be  obtained  from  Eq.  (46). 


17 


SECTION  III 


RESULTS  AND  DISCUSSION 

Because  of  the  complications  in  analysis  of  the  many  coupled  modes, 
only  one-mode  approximation  is  used  in  the  formulation.  The  assumption 
for  fundamental  mode  predominacy  is  admittedly  overly  simplified;  the 
conditions  under  which  this  is  a  valid  approximation  remain  to  be  inves¬ 
tigated.  However,  a  simple  model  sometimes  helps  to  give  basic  under¬ 
standing  of  the  problem. 

Using  the  present  formulation,  response  of  nonlinear  rectangular 
panels  with  all  edges  simply  supported  and  all  edges  clamped  subjected  to 
broadband  random  acoustic  excitation  are  studied.  Both  immovable  and 
movable  inplane  edges  are  considered.  In  the  results  presented,  the  spec¬ 
tral  density  function  of  the  excitation  pressure  S(f)  is  considered  flat 
within  a  certain  region  near  the  equivalent  linear  frequency  f  and  a 
value  of  Poisson's  ratio  of  0.3  is  used  in  all  computations,  unless  other¬ 
wise  mentioned.  Mean-square  modal  amplitudes  and  mean-square  nondimensional 
stresses  for  panels  of  various  aspect  ratios  and  damping  ratios  are  deter¬ 
mined  and  presented  in  graphical  form.  These  graphs  can  be  used  as  guides 
for  preliminary  design  of  aircraft  panels.  The  maximum  mean-square  deflec¬ 
tion  can  be  reasonably  obtained  from  these  figures;  however,  multiple-modes 
had  to  be  considered  for  accurate  determination  of  mean-square  stresses. 

This  has-been  demonstrated  by  Seide  in  Reference  15  for  a  simple  beam 
subjected  to  uniform  pressure  excitation  and  in  Reference  16  for  large 
deflections  of  prestressed  simply  supported  rectangular  plates  under 
static  uniform  pressure. 

Comparison  with  experiment  is  also  given.  It  demonstrated  that  the 
present  formulation  given  remarkable  improvement  in  predicating  RMS 

responses  as  compared  with  using  the  linear  theory. 
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1.  Simply  Supported  Panels 

Figure  2  shows  the  maximum  mean-square  nondimensional  deflection 
versus  nondimensional  spectral  density  parameter  of  acoustic  pressure 
excitation  for  rectangular  panels  of  aspect  ratios  r  =  1,  2,  and  4,  and 
a  damping  ratio  0.02.  It  is  clear  from  the  figure  that  an  increase 
of  r  will  cause  "closing"  of  the  curve.  This  occurs  because  as  r  increases 
the  panel  becomes  less  stiff,  and  the  mean-square  deflection  has  to  be 
finite.  It  can  also  be  seen  from  the  figure  that  the  mean-square  deflec¬ 
tion  of  the  movable  inplane  edges  case  is  approximately  twice  as  that  of 
the  immovable  edges. 

Figures  3  and  4  show  the  maximum  mean-square  bending  stress  and  the 
mean-square  membrane  stress,  respectively,  for  a  simply  supported  square 
panel  with  5  =  0.02.  The  maximum  mean-square  bending  stress  of  the  mov¬ 

able  edges  case  is  approximately  twice  as  much  as  that  of  the  immovable 
edges,  whereas  the  mean-square  membrane  stress  of  the  movable  edges  case 
is  found  to  be  somewhat  less  than  that  of  the  immovable  edges.  In  Figure 
5,  the  maximum  mean-square  nondimensional  stress  (bending  plus  membrane 
stress,  at  the  center  of  the  panel  and  in  the  y-direction)  is  given  as  a 
function  of  excitation  spectral  density  parameter  for  simply  supported  rec¬ 
tangular  panels  of  various  aspect  ratios  and  a  damping  factor  0.02.  Results 
showed  that  the  difference  of  maximum  mean-square  stresses  between  immov¬ 
able  and  movable  edges  is  small  as  compared  with  the  difference  of  mean 
square  deflections  between  the  two  edge  conditions. 

Figure  6  shows  the  mean-square  deflection  versus  forcing  spectral 

density  parameter  for  simply  supported  square  panels  of  different  damping 

ratios.  The  corresponding  maximum  mean-spare  stress  (bending  plus  membrane 

stress,  at  the  center  of  panel),  is  shown  in  Figure  7.  As  it  can  be  seen 

from  the  figures  that  the  precise  determination  of  damping  ratio  from 
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experiment  is  important,  e.g.  stress  increases  by  25-30%  as  t,  is  decreased 
from  0.015  to  0.01  (for  between  5000  to  20000).  Again,  the  difference 

of  maximum  mean-square  stresses  between  inmovable  and  movable  edges  is  small 
as  compared  with  the  difference  of  mean-square  deflections. 

Plots  of  the  equivalent  linear  or  nonlinear  frequency  parameter 

2 

A  versus  mean-square  modal  amplitude  for  simply  supported  rectangular 

panels  of  aspect  ratios  r  =  1,  2,  and  4  are  shown  in  Figure  8.  The 
2 

lowest  value  of  A  corresponds  to  the  linear  case. 

2.  Clamped  Panels 

In  Figure  9,  the  mean-square  deflection  is  given  as  a  function  of 
excitation  spectral  density  parameter  for  rectangular  panels  of  aspect 
ratios  r  =  1,  2,  and  4  and  a  damping  ratio  0.02.  The  maximum  mean- 
square  deflection  of  the  clamped  panels  is  somewhat  much  less  than  that 
of  the  simply  supported.  The  corresponding  maximum  mean-square  nondi- 
mensional  stress  (bending  plus  membrane  stress,  in  the  y-direction  and 
at  the  center  of  the  long  edge)  versus  spectral  density  parameter  is 
shown  in  Figure  10. 

Figure  11  shows  the  mean-square  modal  amplitude  versus  spectral 
density  parameter  of  excitation  for  a  square  panel  of  different  damping 
ratios.  In  Figure  12,  the  equivalent  linear  frequency  parameter  is  given 
as  a  function  of  mean-square  deflection  for  clamped  rectangular  panels 
of  aspect  ratios  r  -  1,  2,  and  4. 

3.  Comparison  with  Experimental  Results 

The  experimental  measurements  on  skin-stringer  panels  exposed  to 
random  pressure  loads  reported  in  References  3  and  4  are  used  to  demon- 
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strate  the  improvement  in  predicting  panel  responses  by  using  the  pre¬ 
sent  formulation.  The  structure  was  a  skin-stringer,  3-bay  panel  as 
shown  in  Figure  13.  The  panels  were  constructed  of  7075-T6  aluminum 
alloy.  Details  of  the  test  facility,  noise  sources,  test  fixture,  and 
test  nssulltLc  are  given  in  Reference  3.  The  important  properties  of 
the  panel  are: 


Length 

Width  between  the  rivet  lines 

Thickness 

Damping  ratio 

Poisson's  ratio 

Young's  modulus 

Weight  density 


a  =  27  inches 
b  =  6.63  inches 
h  =  0.032  inches 
c  *  0.0227 
v  *  0.33 

E  -  9.6  x  106  psi 
p  =  0.1  lb/in3 


The  tests  were  conducted  with  an  overall  sound  pressure  level  (SPL) 
of  157  decibels  (dB)  with  a  range  of  ±  1.5  dB  which  corresponds  to 
an  average  spectrum  level  of  125.26  dB  (see  Table  IV  of  Reference  3 
or  Table  8  of  Reference  17). 

The  central  bay  of  the  3-bay  test  panels  is  simulated  by  a  flat 
rectangular  plate.  The  linear  frequencies  for  both  simply  supported 
{ £q.  (14))  and  clamped  (Eq.  (20))  support  conditions  are  calculated 
and  shown  in  Table  1. 

Test  measurements  and  finite  element  solution  are  also  given 
for  comparison.  Table  1  also  shows  the  equivalent  linear  or  nonlinear 
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Table  1  Frequency  Comparison 


Natural 

Equivalent  linear 

frequency  f 

frequency  fJ57 

Simply  supported  -  Immovable  edges 

71 

321 

-  Movable  edges 

71 

240 

Clamped  -  Immovable  edges 

159 

311 

-  Movable  edges 

159 

264 

Finite  element  (Ref.  4) 

155 

- 

Experiment  (Ref.  3)  126, 


frequencies  at  overall  SPL  157  dB.  Frequency  at  high  intensity  noise 
level  was  not  reported  in  Reference  3.  From  the  results  shown  in  Table  1, 
it  is  clear  that  the  central  bay  of  the  test  panels  did  not  respond  to 
the  acoustic  excitation  as  though  it  were  fully  clamped  on  all  four  edges. 
This  was  also  demonstrated  in  Figures  12  and  17  of  Reference  3  in  the 
sense  that  the  highest  measured  RMS  strains  did  not  occur  at  the  center 
of  the  long  edges.  The  central  bay  of  the  test  panels  actually  behaved 
somewhere  between  fully  simply  supported  and  fully,  clamped  support 
conditions. 

The  acoustic  pressure  spectral  density  S(f)  is  related  to  the 
spectrum  level  L  as 

.  8.41  x  10(L/1°  '  18)  (psi)2/Hz 


S(f)  = 


4  x  10  <L'10  -  8> 


( dynes/cm2 )2/Hz 


(47) 
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A  white  noise  pressure  loading  with  spectral  density  of  S(f)  =  2.824  x  10 
(psi)  /Hz  (or  nondimens ional  spectral  density  parameter  S^,  =  5100),  which 

corresponds  to  an  average  spectrum  level  L  =  125.26  dB,  is  used  in  the  com¬ 
putations.  The  RMS  stresses  (Eq.  (46))  at  the  center  of  the  long  edges  for 
simply  supported  (Eq.  (41))  and  clamped  (Eq.  (43))  boundary  conditions  are 
calculated  and  given  in  Table  2. 


Table  2  Stress  Comparison 
RMS  Stress  (kpsi) 


Linear 

Nonl inear 

Linear 

Nonl inear 

theory 

theory 

theory 

theory 

Simply  Supported 

0.0 

0.58  (Immovable) 

0.0 

3.28  (Immovable) 

0.17  (Movable) 

2.74  (Movable) 

Clamped 

2.17 

1.12  (Immovable) 

6.57 

3.83  (Immovable) 

1.32  (Movable 

4.24  (Movable) 

Finite  element 

2.4 

7.7 

(Ref.  4) 

Measured  (Refs.  3,4) 

on  panel  A 

0.63 

2.2 

B 

0.94 

2.9 

C 

0.78 

2.5 

D 

1.1 

— 

E 

0.84 

2.2 

Measured  average 

0.87 

’2.5 

Table  3  shows  the  RMS  deflections  using  the  present  formulation.  The 
measured  and  finite  element  RMS  stresses  and  RMS  deflections  in  Reference 
4  are  also  given  in  the  tables  for  comparison.  It  demonstrates  that  a 
better  correlation  between  theory  and  experiment  can  be  achieved  when 
large  deflection  geometrical  nonlinearity  effect  is  Included  In  the 
formulation. 
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Table  3.  Deflection  Comparison 

j  <wh)S 


Linear  theory 

Nonlinear  theory 

Simply  Supported  8.0 

1.8  (Immovable) 
2.4  (Movable) 

Clamped  2.7 

1.4  (Immovable) 
1.6  (Movable) 

Finite  Element  (Ref.  4)  3.1 

Measured  (Ref.  4) 

2.0 
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SECTION  IV 


CONCLUSIONS 

An  analytical  method  for  predicating  response  of  rectangular  nonlinear 
structural  panels  subjected  to  broadband  random  acoustic  excitation  is  pre¬ 
sented.  The  formulation  is  based  on  the  Karman -Herrmann  large  deflection 
plate  equations,  a  single-mode  Galerkin  approximation,  the  equivalent 
linearization  method,  and  an  iterative  procedure.  Both  simply  supported 
and  clamped  support  conditions  with  immovable  or  movable  inplane  edges 
are  considered.  Panel  mean-square  deflection,  maximum  mean-square  stress, 
and  equivalent  linear  frequency  at  given  excitation  pressure  spectral 
density  can  be  determined,  and  they  are  presented  in  graphical  form. 

These  graphs  can  be  used  as  guides  for  preliminary  design  of  aircraft 
panels  under  high  noise  environment.  Results  obtained  agree  well  with 
the  experiment.  It  is  suggested  that  further  research  be  carried  out 
with  special  attention  to  employ  multiple-modes  in  the  formulation, 
and  additional  test  data  on  simple  panels  are  needed  for  an  adequate 
quantitative  comparison. 
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Figure  2.  Mean^square  deflection  versus  spectral  density  parameter 
of  excitation  for  simply  supported  panel s,t  =  0.02. 


Figure  3.  Maximum  mean^square  bending  stress  versus  spectral  density  parameter 
of  excitation  for  a  simply  supported  square  panel,  t,  -  0.02 


membrane  stress  versus  spectral  density  parameter  of  excitation 
’  supported  square  panel,  c  =  0.02. 


simply  supported 


gure  7.  Effects  of  damping  on  maximum  mean-square  stress 
square  panel 


Fi9“re  *'  for^lamped  JEST-TS!  SPKtra''  Paramet!n-  °f 


Ill 


Effects  of  damping  on  mean-square  deflection  for  a  clamped  square  panel 


I  I  I  I  I  I 

0.5  1.0  1.5 _ 2J)  2.5  3.0  3.5 

(iT-f 

12.  Frequency  parameter  versus  mean-square  deflection  for 
clamped  panels 
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APPENDIX 


LARGE  AMPLITUDE  RESPONSE  OF  AIRCRAFT 
RANDOM  ACOUSTIC  EXCITATIONS 

SIMPLY  SUPPORTED  ALONG 
RECTANGULAR  PANEL  -  SIMPLY  SUPPORTED  ALONG 

SIMPLY  SUPPORTED  ALONG 
COMMON  AK,AR2,PI,PI2,UN,£NU 
DIMENSION  SPSU(99> 

MOVABL=FL AG  FOR  INPLANE  EDGE  CONDITIONS 
SET  MOVABL=  0,  FOR  IMMOVABLE  EDGES 
SET  MOVA8L=+l,  FOR  MOVABLE  EDGES 
MOV  ABL=0  OR  ♦  !  ? 


FANELS  SUBJECTED  TO  BROADBAND 


ALL  FOUR 
ALL  FOUR 
ALL  FOUR 


EDGES 

EDGES 

EDGES 


MOV  A3L  =  1 

ASGNG=AN  ASSIGNED  SMALL  CONSTANT  FOR  CONVERGENCE  TEST  OF  RMS 
MAXIMUM  DEFLECTION 
ASGNQ= 1 • OE-3 


PI  =  3. 14159265 
PI2=PI*PI 
PI4=PI 2*^1 2 

ZET A  =  D AMP ING  FACTOR  OF  PANEL=C/(C  CRITICAL) 

ZET  A=0 . 020 
UN^POISSON'S  RATIO 
UN=0. 3 

ENU=1.0  -  UN *UN 
9000  rORNAT (12) 

9050  FOR“AT(F10.3) 

NUmSF=MUM8£R  OF  EXCITATION  SPECTRAL  DENSITY 
800C  FJRMAT(43H  HOW  «ANY  EXCITATION  SPECTRAL  DENSITY  ARE  THERE?# 14) 
READ  (5,9000)  NUMSF 
WRITE  (6,8000)  NUMSF 
-RITE  (6,9100) 

DO  100  \’SF=1, NUMSF 

8100  FORMAT ( 23H  WHAT  ARE  THEIR  VALUES?  ) 

HEAD  (5,9050)  S0SD(N3F) 

-RITE  (6,8150)  SPSD(NSF) 

8150  FORVAT(5X,E15.o) 

SPSD=NOND IMENSIONAL  FORCING  EXCITATION  SPECTRAL  DENSITY 
PARAMET£R=Sf(F)/((RH0**2)  *  (H**4>  *  (DFACTOR**1.5)) 

IOC  CONTINUE 

C  NAh=NUMPEP  OF  ASPECT  RATIOS 
8200  FORMAT (24H  HOW  MANY  ASPECT  RATIOS?, 14  ) 

READ  (5,9000)  NAR 
WRITE  (6,8200)  NAR 
DO  900  NTIME=1,NAR 
C  AR= ASPECT  RATIO  OF  PANEL= AL/3L 
C  AL=LENGTH 
C  BL=WIDTH 

READ  (5,9050)  AR 

Af?2=AF*AR 

AR4=AR2*AR2 

C  FLkD02=LAMBD AO  SQUARE  -  NONDIMENSIONAL  LINEAR  FREQUENCY  SQUARE 
FLMD02=P14*(1.0  ♦  AR2)*(1.0  ♦  AR2)/AR4 
FLMPO=  SQRT(FLMT)02) 

C  BET AP=NON LINEARITY  COEFFICIENT 
C  EET AC=NONLINEARITY  COEFFICIENT 

BFTAP=0.75*PI4*ENU*(1.0  ♦  1.0/AR4) 

BETAC=1.5*PI4*(I.O  +  UN*AR2  ♦  AR2*(AR2  ♦  UN))/AR4 
C  FOK  MOVABLE  INPLANE  EDGES,  ^ETA=BETAP 
C  FOR  IMMOVABLE  INPLANE  EDGES,  BETA=3ETAP  ♦  BFTAC 
BET A=BETAP  ♦  BETAC  39 


n  ci  o 


IF  (M0VA3L  .fO.  1)  3£TA=3ETAP 
C 

UO  ^00  NSF=1,NUV3" 

3F=SPSD(N3F) 

WRITE  (6,9100)  AR,  SF,QETA 
9100  FOkwAT(//17HlA3°ECT  R AT  10= AP  =  ,F6 . 2/, 

1  47H  N0NPIMEN3IU*!AL  EXCITATION  SPECTRAL  D£NSITY=SF=, 

2  E13.5/31H  NONLINEAR  I TY  COEFF ICIENT=BETA=,E15.6) 

IP  (MOVABL  .29.  0)  WRITE  (6,8300)  ZETA 

8  300  FOR  FAT (46H  SIMPLY- SUPPORTED  WITH  IMMOVABLE  INPLANE  EDGES, 

1  / 3 4 M  DAMPING  FACTOR=C/C  CR IT ICAL=ZETA=,F7 .4) 

IF  (MOVAaL  .  EQ.  1)  WRITE  (6,8310)  ZETA 
3310  FORMAT(44H  SIMPLY- SUPPORTED  WITH  VOVAPLE  INPLANE  EDGES, 

1  / 34H  DAMPING  FACTOR=C/C  CRITICAL=ZETA=,F7.4) 

C  ITRR=0  CORRESPONDS  TO  THE  LINEAR  SOLUTION 
ITER=0 

WRITE  (6,9150) 

9150  FORMAT (//41H  *** SMALL-DEFLECT  ION  (LINEAR)  SOLUTION*** ) 

C  G02=ME AN  SQUARE  MAXIMUM  DEFLECTION  -  LINEAR  STRUCTURAL  THEORY 
902=32. 0*SF/(',14*Z2TA*FL.MDQ*FLMDQ2) 

C  QO=RCOT  MEAN  SQUARE  MAXIMUM  DEFLECTION 
«0=  SQRT(OQ2). 

WRITE  (6,9200)  FLMD02,FLMD0, Q02, QO 
9200  FOR  MAT ( 4X, 3 1 H  FRE9.  PARAMETER  SQUAR£=FLMD02=,E15. 7, 

1  5X,23H  FREQ.  PA,’AWETEP=FLMDO=,E15.7/, 

2  8x, 27H  MEAN-SQUArE  AMPLITUDE=902=,E15.7, 

3  10x,18H  RMS  AMPLITUDE=QO=,F7.3) 

C  COMPUTE  R VS  STRESSES  OR  STRAINS  EASED  ON  LINEAR  STRUCTURAL  THEORY 
aPITE  (6,9250) 

9250  FORMAT (/3X,9H  LOCATION, 20X, 19H  MEAN-SQUARE  STRESS, 19X, 

1  I1H  RMS  STRESS/ 29X, 9H  X-C3M.PS 6X, 9H  Y-COMPS. , 14X, 

2  9P  X- COMPS. ,6X,9H  Y-COMPS.) 

CALL  STRESS(0.0,0.0,Q02,MDV ABL, ITER) 

CALL  STKESS(0.5,0.0,Q02,KOV ABL, ITER) 

CALL  STRESS(0.0,0.3,Q02,MOVABL,ITER) 

CALL  STRESS(0.b,0.5,Q02,MDVABL, ITER) 

C  STORE  Q02,  QO  AND  FLdDO  FOR  CONVERGENCE  TEST  LATTER 
Q2P V=OQ2 
QPV=QO 

FLMDPV  =FLMDO 

START  UN  ITERAHON  -  ITER=ITERATION  COUNTER 
I TER=1 

10  CONTINUE 

FLrtD2=FLMD02  +  3.0*8ETA*Q2PV 
Q2=32.0*SF/(PI4*ZErA*FLMDO*FLMD2) 

FLMD=  SQR?( FLVD2 ) 

9=  S3RT(q2) 

C  CHECK  FOR  CONVERGENCE 

DEF=(FLMD-FLMDrV)/FLMC 

DEG=(Q-QPV)/Q 

DFQ=  ABS(DEO) 

C  COfVERGENCE  TEST,  IF  SATISFIED  go  TO  12  FOR  STRESS  COMPUTATION 
IF  (OEq  .LE.  ASCNQ)  GO  TO  12 
C  STORE  FLMD,  Q2  AND  0  FJa  CONVERGENCE  TEST 
FLMDPV=FLMl> 

92PV=  32 

C  SAWE  COMPUTATION  TI^E  AFTER  50  AND  100  ITERATIONS 
IF  (ITER  .EO.  50)  Q2PV  =  0 .5^  Q*9  ♦  QPV*QPV) 


OCiO 


IF  (ITSP  .c.n.100)  Q2FV=0.5*(Q*Q  ♦  QPV*QPV> 
gPV=0 

1TER= ITFR  ♦  1 
GO  TO  10 
12  CONTINUE 

WRITE  (6,9400)  ITER 

940G  2 Ok  FAT ( // 44  u  ***LARC2-DE1-’LECT  ION  (  NONLINEAR)  SOLUTION***/, 

1  -4X,  loH  CONVERGENCE  AT,14,13H-TH  ITERATION) 

Mil  TE  (6,94 SO)  FLMD2,  FL.MD,  Q2,  0 
9450  r'0kPAT(4X,3CH  FRE7.  PARAMETER  SQUARE  =  FD*D2=,E15.7, 

1  bY,22H  PREJ.  PAPAM2TSR=FLMD=, E15.7/, 

2  8x,  26H  VEAN-SOUAkE  AMPLITUDE=72=,S14.6, 

3  llx,17S  RMS  APPLI  rL'DE=Q  =  ,F7.  j) 

WRITE  (6,9500)  DEP,D2'2 

950C  fOkf,AT(?OX,4Hp*EF=,sl5.6,bX,4HDEO=,E15«6) 

C  CO!  PUTE  STRESSES  A:'  LJRGE  AMPLITUDE 
*RITS  (6,9250) 

CALL  3TR2S3(0.0,0.0,Q2,VUVA?L,ITER) 

CALL  STRES3(0.6,0.0,G2,MOVA3L, ITER) 

CALL  STPE33(0.0, G.b,G2, MOV ABL, ITER) 

CALL  STRESS(0.5,G.S,92,V0VAPL,1TER) 

IP  ( 4  .02.  c.O)  G3  TrJ  700 
500  CONTI "UK 
700  COMIN"E 
9  0  C  CfK.TINUF 

3T0'1 
2  ML 

3" 3 ROUTINE  ?T? 2  3S(  X  A,  VE ,rj2,  MQVAPL,  ITER  ) 

C'JMM.jn  AH,  A°2,?  I,*’I2,Uii,£NU 

C  THIS  IUSkHUTIVL-  C1MPITS  3  THE  RMS  STRESSES  IN  THE  X  (3X)  AND  Y 
C  (5V)  DIRECTIONS  -  rHs.SE  STRESSES  APE  .'JONDIhENSlONAL 
C  XA-X/A 
C  Y" - Y/ A 

C  Oi-’  EAu  SODA  RE  "2  '1 
CX  =  OjS ( ?i*X  A ) 

CY=CJS(P1*Y?) 

J2X=COS(i.O*PT*XA) 

C2Y=Cns(2.0*oI*YE) 

7X  =  0.5*?I‘>*(ON  *  l.D/Art2)*CX*CY/2»lU 

31X=0.125*P!2*C2Y/AR2 

i)  2  a  =  0  .  1 2 5 * ?  1 2 * ( fi Y  +  I.0/AS2  )/EN’H 

IF  (.VO VAHL  .2°.  1)  ^4,X  =  0.0 

D v=0. 5  *°12* ( 1 . 0  i-  UN/  AP2)*CX*CY/SNU 

:;lY=o.  i?b*r  T2*c?x 

J2  i=C.12j>*P72*(U:I  >  A R2  )/(  A?2 *ENU ) 

I  p  (MOVAEL  .£•).  I)  0  2  v=  0.  0 
DXX-OIX  ♦  02 X 

c  irf-r,  =  n  CORRESPONDS  TG  r»F.  SMALL  DEFLEC  T  ION  LINEAR  THEORY 
IF  (ITER  .00.  0)  ">XX=0.O 
SXi-  R=PX*9X*02 
SXM2=3.0*rXX*CXX*92*  i? 
jX2=3Xb2  +•  SX'-'s 
Ovv=')iY  +  02Y 
IP  (ITER  .EC.  0)  0  Y v  =  0 . 0 
J  Y  2  2  =  P  Y  *  J  v  *  0  2 
SY S 2=3.n*PYY*rYY*Q2*G2 
■>v2-3VF2  *  > Y**2 

SX,  SY=''I‘W.N.CICVLFSS  root  I  FAN  SO’UREH  STRESSES  ( STRESS*! **2 ) 

/  (  E  *  :i  *  *  2  ) 

F  = '  L='V I DT F  OF  •’AflEL 
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c  E=FODULUS  OF  ELASTICITY 
C  H=TH1CKNESS  OF  PANEL 
SX=  SQRT(SX2) 

SY=  SQRT{SY2) 

WRITE  (6/9000)  XA/SX2/SY2/SX/SY 

9000  FORFAT{3x/5h  X/A=/ F5 . 7, 5X/ 6H  TOTAL/3X/2E15.6/ 8X/2E15.6) 
9200  FORMAT ( 3x/ 5h  Y/A=/F5 . 2/ 5X/ 8H  9ENDING/1X/2E15.6/ 8X/2E15.6) 

SXB  =  SQRT (5X32) 

SXK=  SQRT(SXK2) 

SYB=  SQRT(SY82) 

SYM=  SQRT(SYM2) 

WRITE  (6/9200)  YB/SXB2/ SYB2/SXB/SYB 
WRITE  (6/9400)  SXM2/SYM2/SXM,SYM 
9400  FORMAT (18X/9H  MEMBRANE/2E15.6/8X/2E15.6) 

RETURN 

END 
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o  o  o  o  on  no 


C  LA :•■<?£  AJ*t*LITUnii  RESPONSE  Cl  AIRCRAFT  PANELS  SUBJECTED  r  '1RC„">[  A!.n 
RANDOM  ACOUSTIC  EXCITATIONS 

CLAPPED  ALONG  ALL  FOUR  EDCc.3 
k  EC  TAN  SUL  i  R  PANEL  -  CLAVPZD  ALONG  ALL  FOUR  E0G&3 

CLAMPED  ALONG  ALL  FOUR  EDGES 

COM  VON  A.i2/oI/?I2/UN/i.*,J/°H10/PH01/PE20/PHll/PI:02/F!,4il/oK12 
DIMENSION  3?3D(99) 

F  OVABL-FLAG  rC.<  IM-LANE  i'lGi  CONDITIONS 
SET  MOV  AST.  =  0,  FOR  I M  VOVA  OLE  EDGES 
3 ET  M  0  V  A  3  L  =  ♦  l ,  FOR  MOVABLE  EDGES 
MO  V  A  EL =0  CR  >1  ? 

MOV  APL  =  0 

C  ASGNQ=AK  ASSIGNED  SMALL  CONSTANT  FOR  CONVERGENCE  TEST  OF  RES 
C  MAXIMUM  DEFLECTION 
ASGN  1-1.0E-3 
C 

Pl=3. 14159265 
P  1 2  =P  I  *P  I 
PI4=?I2*PI2 

C  ZETA=r'AFI-  ING  FACTOR  OF  ?AN5L=C/(C  CRITICAL) 

ZET  A  =  0 . 020 

C  UN-POISSON'S  RATIO 
UN-0.3 

E  NU  =  1 . 0  -  U  N * U  N 
900C  FORMAT (  12  ) 

905*'  FORMAT  ( F 1 0 . 2  ) 

C  NUN  SF  =  , MUMPER  OF  EXCITATION  SPECTRAL  DENSITY 
9000  rORKAT(43H  MOW  *kVl  EXCITATION  SPECTRAL  DENSITY  ARE  T^ERS?,I«) 
REaD  (3,9000)  NU.V3F 
-RITE  (6,3000)  M U M $ E 
..RITE  (6,9100) 

DO  100  NSF  =  1 , NJVS7 

3100  FORMAT ( 23H  WHAT  ARE  THE  I R  VALUES?  ) 

'READ  (  3,  TO 50  )  SPSD(  MSF) 

WRITE  (6,3150)  SPSD(NSF) 

3150  rORMAT(5X,E15.6) 

C  SPSD=NOM DIMENSIONAL  FORGING  EXCITATION  SPECTRAL  DENSITY 
C  PARAMETER=SF(F)/((RH3**2)  *  (H*M)  *  (OF ACTOR**l . 5  )  ) 

10G  CONTINUE 

C  M Ak  =NUMBER  OF  ASPECT  RATIOS 
8200  FQRMAT( 24H  HOW  *'AN Y  ASPECT  RATI0S?,I4  ) 

READ  (5,9000)  NAS 
WRITE  (6,8200)  NAR 
DO  900  NTIME=1,NAH 
C  AR=ASPECT  RATIO  OF  PANEL= AL/EL 
C  AL=LENGTH 
C  BL-WIDTH 

READ  (5,9050)  AR 

AR2=AR*AR 

AR4=AR2*AR2 

C  FLKD02=LA  KBQ  AQ  SQUARE  -  NONDI KENS  I CNAL  LINEAR  FREQUENCY  SQUARE 
FLMDO2=16.0*PI4*(3.0  <■  2.0*AR2  ♦  3.0*AR4)/(9.0* AR4) 

FLMDO=  SQRT (FLMD02  ) 

C  BETAP=NONLINEARITY  COEFFICIENT 
C  BET AC=NONLINE ARITY  COEFFICIENT 
PH10=1.0 
PHO 1=1 .0/ AR4 
PH20=0 .0625 

PHI 1=2 . 0/ ( 1. 0  ♦  2. 0*AR2  +  AR4) 

PH02=0.0625/AR4  43 


O  o  o 


Pf’21  =  1.0/(  16. 0  *  P.O*AK?  +  A.'<4) 
i-'*12=1.0/(l-0  *  3.0*AP2  +  15.0*AH4) 

PHSUM=P?!10  *  P  iO  1  «•  P'I20  *  P.m  *  P M 0 2  ♦  0.r*(P'<21  ♦  PP.12) 
.  i  E  T  A  P  =  1 . 0  *  P  I  4  *  Z  N ! )  *  P  M  i : j  M  /  J  .  0 

OFT  AC  =  1 .  j*P  I  4* ( 1 . 9  ♦  U;.*AR2  ♦  APz*(Ak2  +  P:. ) )  /  A  ^  < 

C  FT.  MOVALLE  INPLmWK  iiDG23,  rrTA=PFTA? 

C  FOi  I  M  M  0  V  A  a  L  v  lKPLANr.  GOGHS,  :3ZTA  =  EET  AP  +  B  H  7  A  C 
SHI  A  =  BET  AP  ♦  °ZT  AC 
lc  (MOVASL  .  HQ.  1)  q " T A  =  r F T A P 

c 

DO  50P  NSF  =  1,:!'JV3F 
i  F  -  5  P  S  D  (  :J  S  F  ) 

XPITE  (  6,  9 1  0  J  )  Ak,3F,3GTA 
910"  FDRVATC//17M1  ASPECT  E  AT  TO  =  AD=,F  4  .  2/, 

1  47  K  NONIU'-'ZKSIONAL  HXCITATi'Ji.  SPECTRAL  DEN  J  IT  Y=  "■? =  , 

7  „13.5/31’J  NONLINEARITY  C  0 H  t  F I  C I H ;<  T  =  °  E T  A  =,  E 1 5 . 6  ) 

if  (mOVAJL  .so.  (j)  ..’RITE  (o,8300)  ZEXA 
9300  F 1-1  A  T  (  4  7!!  Ci-A  -lopO- GU  PPlETET  V I  ^  H  IMMOVABLE  INFLAME  EDGES, 
i  /34»  0  A '  V I N  C  FACTOR=C/C  CRITIC AL  =  ZETA  =  /F7. 4) 

IF  (MOV A  PL  .El.  1)  WHITS  (o,8310> 

9310  rMr'-'ATdj'!  C^  5  M  PHD- G'JPP  3«TED  WITH  MO  V  A  CLE  I  ‘.'PLANE  EDGES, 

1  /  3  4  H  "ANVlNC  FACT  0  R  -  C  /  C  CR  IT  I C  AL  =Zc.T  A  =  ,  F7 . 4) 

c  it:.k=o  copkhs°jnds  tj  tue  li:*ea°  solution 

I  T  E  R  =  C 


WRITE  (5,9150) 

9150  PQRVATC//4I"  ***3V«LL-D EFFECT  ION  (LINEAR)  SOLUTION***) 

C  CO:  =  fcSAN  SQUARE  N  AX  I  MUR  REFLECTION  -  LINEAR  STRUCTURAL  THEORY 
C02=3  2 . 0  *Sr /( 3 1 . 0*ZET A*FLMDO*FLMDQ2 ) 

C  GO-ROOT  MEAN  SQUARE  MAXIMUM  DEFLECTION 


JO  =  SORT ( QC2 ) 

WRITE  (6,9200)  ^L'-'DOW/FLMD  J/UOi,  JO 


920a 

1 

A 

2 

3 

C  CJi-.PU 
9250 

1 

2 


F0RMAT(4:</J1.I  FPHQ.  PARAMETER  SOU  AkE=FLMD02=, F 15 • 7, 

3X, 230  FREQ.  ?APAMETE9=FLMDO=,El 5 .7/, 

8x,27"  “EA., -SQUARE  AVPLI TU'JE  =  Q02  =  ,E15. 7, 

1 0  x  / 1  9  H  RMS  AMi,LIT'Jr'S  =  QO  =  /F7.3) 

TE  RMS  STRESSES  OR  STRAINS  BASED  ON  LINEAR  STRUCTURAL  THEORY 
wN  I  TiZ  (5,9250) 

FORMAT  (  /3X,9'H  L0CAT  I  ON,  20X,  1  9H  VEAN-SQUARE  STRESS,  PX, 

1 1 H  RVS  3TRE3S/29X,  9H  X-C0^?S.,6X,9H  Y-CJ  VPS.  ,  1 -1 X, 

9H  Y-CiJMPS.,5X,9;I  V-COMPS.) 

CALL  3 TRESS (0.0,0.0,302, MOV AHL, ITER) 


CALL  STRESS(0.5,0.0,  J02,  MOVAF-L,  I  TER) 

CALL  S T R£SS(0.0, 0.5,00 2, MOVAoL, ITER) 

CALL  ST  PE  SS(0. 3,0. 5,  C  02,  MOVASL,  HER) 

C  STORE  w02,  QO  AMD  FLMDO  FOP  CONVERGENCE  TEST  LATTER 


Q2P V=Q 02 


QPV-QC 


FLMDPV  =FLMDO 


START  OK  ITERATION  -  I  TER  =  ITER AT  I  ON  CHUNTER 


ITE°=1 

10  CONTINUE 

r LKD2=rLMD02  ♦  3. 0*°ETA*Q2PV 
J2=32.0*SF/(81.O*ZETA*FLMDO*FLMD2) 
FLMD=  SQrT (FLVD2) 

Q=  SQPT(Q2) 

C  C°:.rK  FOR  CONVERGENCE 

DEF=(FLMD-FLMT)PV)/FLMD 
OEQ=(Q-QPV)/Q  44 


DEQ  =  ABS(DKQ) 

C  CONVERGENCE  TEST,  IF  SATISFIED  GO  TO  12  FOR  STRESS  COMPUTATION 
IF  (DEQ  .LE.  ASGNQ)  GO  TU  12 
C  STORE  FLMD,  Q2  AND  Q  FOR  CONV ERGENCF.  TEST 
f LNUPV=FLMD 
Q2?V=Q2 

C  SAVE  COMPUTATION  TI^E  AFTER  50  AND  100  ITERATIONS 
IF  (ITER  .£0.  aO)  Q 2F V=0. 5* ( G*Q  ♦  QPV*GPV) 

IF  (ITER  .ED. 100)  Q2PV=0. 5* ( QPV*GPV) 

QP  V  =Q 

iT£R=ITER  +  1 
GO  TO  10 
12  CONTINUE 

WRITE  (6,9400)  ITER 

9400  FORPAT(//44H  *** LARGE-DEFLECT  ION  ( NONLINEAR)  SGLUTIOw***/ , 

1  4X, 15H  CONVERGENCE  AT, 14, 1 3H-TH  ITERATION) 

-RITE  (6,9450)  FLMD2,FLMD,Q2,Q 
9450  FORMAT(4X,30H  FREQ.  PARAMETER  5QUARE=FLM02=,E15.  7, 

1  5X,22P  FREQ.  PASAMETER=FLWD=,£15.7/, 

2  8x,  26H  MEAN-SQUARE  A/,PL  ITUDE=Q2  =  ,  El  4 .6, 

3  1 lx, 17H  RMS  AMPLITUOS=Q=,F7.3) 

»RIT£  (6, 95C0)  DEF, DEO 

9500  for VAT( 20X,4HDEF=, El  5. 6, 5X,4KDEQ=,E15. 6) 

C  COMPUTE  STRESSES  AT  LARGE  AMPLITUDE 
WRITE  ( 6,9250 ) 

CALL  3TRE3S(0. 0,0.0, Q2, uOV  A PL, I  TER) 

CALL  STRESS ( 0. 5, 0. 0, G2, MJVABL, ITER) 

’ALL  S  TRESS  (  0. 0,  0. 5,  G  2,  MOV  A  3L,  ITER) 

CALL  STRESS(0.5,0.5,r2, uOV A°L,iTSR) 

IF  (Q  .GE.  5.0)  CO  TO  700 
500  CONTINUE 
700  CONTINUE 
900  CONTINUE 
STOP 
END 

SU’t  ROUTINE  STR  c.S3  (  X  A,  YE,Q2,MOVAPL,  I  TER  ) 

COMMON  AR2,  PI,  P 12,  UN,  SN",  PH10,  H»Q  1,  PK2Q,  PHI  1,PH02,P'*21,  PHI  2 
C  THIS  SUBROUTINE  COMPUTES  THE  RMS  STRESSES  IN  TilE  X  (SX)  AND  Y 
C  (Si)  DIRECTIONS  -  THESE  STRESSES  AkE  NONDIMENSION AL 
C  X A=X/ A 
C  Y9=Y/B 

C  Q2=M£ AN  SQUARE  OF  Q 

CX-=CQS(2.0*P1*XA) 

CY=COS(2.0*PI*Y3) 

C2X=CCS(4.0*PI*XA) 

C2Y=COS(4.0*PI*7b) 

UX=0.5*PI2*(  CXM1.0  +  '  «)/AR2  +  UNi*(1.0  +  CX)*CY  )/£NU 
FAC1=0.125*PI2*AR2 

C1X=F4C1*(PH01*CY  +  P  HI 1 *CX*C  Y  +  4.0*PH02*C2Y  + 

1  PP2 1*C2X*CY  ♦  4.0*PF12-*CX*C2Y) 

D2X=3.0*PI2*(1.0  ♦  UE*AR2)/(32.0*AR2*ENn) 

IF  (MQVA3L  .EG.  1)  D2X=0.0 

DY=0.5*PI2*(  (l.O+CX)*CY  ♦  US*CX* ( 1 . 0+C Y )/ AR2  )/ENU 
r AC2-0.!25*PI2 

D1  Y=FAC2*(P',10«'CX  +  4 . 0 *PH2 0* C2 X  +  PH11*CX*CY  ♦ 

1  4.0*PH21*C2X*CY  ♦  PM12*CX*C?Y) 

D2Y  =  3.0*PI2*(  AR2  *  U;i)/(32.0*AR2*ENU) 

IF  (MDVABL  . EQ.  1)  D2Y=0. 0 
DXX-D1X  +  D?X 

C  ITER=0  CORRESPONDS  To  THE  SMALL  ^FLECTION  LINEAR  THEORY 


IF  (ITER  .EQ.  0)  DXX=0. 0 
SXB  2=DX*DX*Q2 
SXM2=3.0*DXX*DXX*Q2*Q2 
SX2=SXB2  ♦  SXM2 
DYY=D1Y  +  D2Y 
IF  (ITER  .EQ.  0)  0 YY=0. 0 
SYB2=DY*DY*Q2 
SYM2=3.0*DYY*DYY*Q2*Q2 
SY2=SYB2  ♦  SYM2 

C  SX,  SY=DIMENSIONLESS  ROOT  KEAN  SQUARED  STRESSES  ( STRESS*B**2) 

C  /(E*H**2) 

C  B=bL=WIDTH  OF  PANEL 
C  E=KODULUS  OF  ELASTICITY 
C  H=THICKNESS  OF  PANEL 
SX=  SQRT(SX2) 

SY=  SQRT(SY2) 

WRITE  (6,9000)  XA,SX2,SY2,SX,SY 

9000  FORMAT ( 3x, 5h  X/ A=,F5. 2,5X>6H  TOTAL,3x,2E15.6,8X,2E15.6) 

9200  FORMAT ( 3x, 5h  Y/A=,F5. 2,5X,8H  BENDING, lx, 2E15. 6, 8X, 2E 15. 6) 

SXB=  SQRT(SXB2) 

SXM=  SQRT (SXM2 ) 

SYB=  SQRT (SYB2) 

SYM=  SQRT(SYM2) 

WRITE  (6,9200)  YB,SXB2, SYB2,SXB,SYB 
WRITE  (6,9400)  SXH2,SYM2,SXM,SYM 
9400  FORMAT ( 18x,9H  MEMBRANE, 2E15. 6, 8X, 2E15. 6) 

RETURN 

END 
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